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ABSTRACT 

We consistently analyse for the first time the impact of survey depth and spatial resolution 
on the most used morphological parameters for classifying galaxies through non-parametric 
methods: Abraham and Conselice-Bershady concentration indices, Gini, M20 moment of 
light, asymmetry, and smoothness. Three different non-local datasets are used, ALHAMBRA 
and SXDS (examples of deep ground-based surveys), and COSMOS (deep space-based sur¬ 
vey). We used a sample of 3000 local, visually classified galaxies, measuring their morpho¬ 
logical parameters at their real redshifts (z ~ 0). Then we simulated them to match the redshift 
and magnitude distributions of galaxies in the non-local surveys. The comparisons of the two 
sets allow to put constraints on the use of each parameter for morphological classification 
and evaluate the effectiveness of the commonly used morphological diagnostic diagrams. All 
analysed parameters suffer from biases related to spatial resolution and depth, the impact of 
the former being much stronger. When including asymmetry and smoothness in classification 
diagrams, the noise effects must be taken into account carefully, especially for ground-based 
surveys. M20 is significantly affected, changing both the shape and range of its distribution at 
all brightness levels. We suggest that diagnostic diagrams based on 2 - 3 parameters should be 
avoided when classifying galaxies in ground-based surveys, independently of their brightness; 
for COSMOS they should be avoided for galaxies fainter than F814 = 23.0. These results can 
be applied directly to surveys similar to ALHAMBRA, SXDS and COSMOS, and also can 
serve as an upper/ lower limit for shallower/deeper ones. 
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1 INTRODUCTION 


Morphology is one of the main characteristics of galax¬ 
ies, and the morphological classification has been central to 
many advances in the picture of galaxy formation and evolu¬ 
tion. Different correlations between morphology and other galaxy 


lar mass (e.g.. 

Deng 2013h. 

colour (e.g.. 

Strateva et alj 200 ll: 

Hogg et al, 2003 

; Bell et alJl2003l: Baldrv et akl 2004; Weiner et al. 

2003; 

ICirasuolo et alj 20051: iMelboume et al. 

12007 ; ICassata et al. 

2007^; 

Povic et alj 2013t). luminosity (e.g.. 

Blanton et al. 20031: 

Kelm et al. 2005h. environment (e.g.. Cassata et alj 2007h. black 

hole mass (e.g..lKormendv & Richstone 

19951: McLure et al.l2000l: 

Graham et al. 

l2001al; Marconi & Hunt 

2003). nuclear activity 

(e.g.. 

Adams 

Il977: HeckmanI 1978: Ho. Filinnenko. & SargentI 

199i 

Kauffmann et al. 200.31; 

Choi et alJl2009l; IPierce et aklbOOTl: 

Povic et aLll2012l; Gabor et al 

l2009h. a 

nd X-ray properties (e.g.. 

Hickox et alj|2009l: Povic et alj2009allbb 


The methods of morphological classification of galaxies 


can be separated into three moups: i) visual (e.g . 
jNajr^^Al^hai^ 201 u [Baillard et alj 1201 ll : 


200S; 


2011 


Lintott et al 


____ iLintott et aL. 

Willett et al.l l2013li . for classifying the nearby and well 


(e.g. de Vaucouleurd 19481; Sersid 1963 

; Peng et alj 20021 

20Kj; 

ISimard et alJl2002Ll201 lUde Souza et alJ 

20041; iBarden et al.1 

2012) 


for fitting the galax' 

M (e.g. Kellv & McKav 

Sooi 

2005:lNgan et al.1 

I 2 OO 9 I; lAndrae et al. 

I 2 OI 1 Jbl: l.limenez-Teia et al. 

I 2 OI 2 I). thev per- 


form the decomposition of well resolved galaxies and provide the 
properties of different structures, and finally, hi) non-parametric, 
which does not assume any particular analytic model. 

The non-parametric methods are based on measuring the 
different galaxy quantities that correlate with the morphologi- 


ties (e.g. Humasonlll93ll: 
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ods are less time-consuming in comparison with other meth¬ 
ods, and can provide an easy and fast separation between reg¬ 
ular and irregular (or disturbed) sources and/or early- and late- 
type galaxies, down to intermediate redshifts (~ 1.5), or higher 
if dealing with space based, or good seeing ground-based data. 
Over the past years different morphological diagrams have been 
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of the previous works showed that these parameters can de¬ 
pend significantly on e.g. t he aperture definition for measur¬ 
ing the galaxy flux/radius l lStrateva et alj l200ll : I Graham et al.l 


noise ratio (hereafter S/N; Conselice et al, I 2 OOOI; iGraham et al. 
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served in relation to the data quality, spatial resolution and depth. 
Since a reliable morphological classification is essential for galaxy 
formation and evolution studies, it is crucial to disentangle how 


strong is this impacQ on each morphological parameter when deal¬ 
ing with different data sets. In most of previous works this obser¬ 
vational bias was analysed either for a particular parameter, either 
for a particular observational condition (survey), or for a particu¬ 
lar sample of galaxies. In this work we present, for the first time, a 
systematic study, where in a consistent way the analysis was carried 
out for all parameters, on large sample of galaxies, and in relation 
with different observational conditions. 

We analysed the observational bias that might affect each of 
the six commonly used morphological parameters: Abraham con¬ 
centration index, Gini coefficient, Conselice-Bershady concentra¬ 
tion index, M20 moment of light, asymmetry index, and smooth¬ 
ness. We studied how this bias depends on the spatial resolution and 
magnitude/redshift distributions of galaxies, and how it affects the 
diagnostic diagrams used to classify galaxies. We used a visually 
classified sample of local galaxies and simulated them to map the 
observational conditions of three different ground- and space-based 
deep surveys: the Advanced Large Homogeneous Are a Medium 


Band Redshift Astronomical survey ijMoles et al.| 

2008, ALHAM- 

BRA), the Subaru/XMM-Aewfon Deep Survey i 

Furusawa 

2008 

SXDS). and the Cosmos Evolution Survey dScoville et al 

2007 


COSMOS). With this analysis our main goal is to measure a set 
of morphological parameters in the real (local) conditions and in 
the simulated conditions of non-local surveys in order to study how 
the spatial resolution and data depth affect each parameter and the 
commonly used morphological diagrams. We emphasize that is out 
of our aim to classify the galaxies in the three selected non-local 
surveys, neither reclassify local galaxies once they were scaled to 
map the conditions of ALHAMBRA, SXDS, and COSMOS. Com¬ 
paring the parameters of local galaxies measured before and after 
shifting them in redshift and magnitude, we were able to put con¬ 
straints on the main diagrams, observing how the position and the 
shape of the regions typical of early- and late-type galaxies change 
in the local and non-local conditions. Moreover, analysing the mor¬ 
phological diagrams in ground-based and space-based surveys, and 
at different magnitude cuts, we quantified how strong is the impact 
from spatial resolution and survey depth, respectively. The results 
obtained in this work can be applied to surveys similar to ALHAM¬ 
BRA, SXDS, and COSMOS. 

The paper is organized as follows: in Sec.l^we described the 
used local and high-redshift samples and the corresponding sur¬ 
veys. The applied methodology is described in Sec. The results 
obtained for local and simulated galaxies are showed in Sec. |4] 
In Sec. |5] we discussed how the impact from data depth and spa¬ 
tial resolution affects the morphological parameters and commonly 
used diagnostic diagrams for galaxy classification. Finally, we sum¬ 
marised our results in Sec.|^ 

We assume the following cosmological parameters through¬ 
out the paper: Da = 0.7, Dm = 0.3, and Ho = 70 km s“* Mpc“*. Un- 
less otherwise speci fied, all magnitudes are given in the AB system 
( lOke & Gunnlll983h . 


^ Hereafter, when using the expression ’observational bias’ we will refer 
to the impact of spatial resolution and data depth on analysed morphological 
parameters. 
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2 THE DATA 
2.1 Local sample 

To test how the observational bias affects the analysed morpho¬ 
logical parameters, we used a sample of 3000 local galaxies at 
0.01 < z ^ 0.1 (with a mean redshift of 0.04), observed in the Sloan 
Digital Sky Survey (SDSS) Data Release 4 (DR4) down to an 
apparent extinction-corrected magnitude of g<16, and visually 
classified by Nair & Abraham (2010; hereafter N&A), using 
the g and r bands. The number of local galaxies is selected as a 
compromise between the computing time and the classification 
accuracy, since the computing time to train the Support Vector 
Machine (SVM) through the galSVM code used in this work 
(see Sec. [3} is totally dependent and very sensitive to the s ize of 
the training data set dHuertas-Companv et al.l 1200 SL l2009ll . The 
galaxies were selected randomly out of ~ 14,000 sources contained 
in the N&A catalogue, making sure that the selected sub-sample 
is represent ative in terms of the general properties of the whole 
data set (see lPovic et alj2013L and their Fig. 3): g band magnitude, 
redshift, g-r colour, morphological classification, and inclination 
in the case of the late-type galaxies. On the other hand. Figure [T] 
shows the redshift, g-band magnitude, and N&A morphological 
classification (T-Type in their work) distributions of the selected 
local sample. As can be seen from the T-Type histogram, the 
selected sample occupies all range of morphologies, from elliptical 
to irregular galaxies. In this work, all analyses were performed 
dividing galaxies into three morphological groups: 

1) early-type galaxies (hereafter ET) - with T-Type <0, 
including elliptical (cO, EO, E+), lenticular (SO-, SO, S0+) and SO/a 
galaxies from the N&A classification, 

2) early spirals (hereafter LT_et) - with 0< T-Type <4, 
including Sa, Sab, Sb, and Sbc, and 

3) late spirals and Irregular galaxies (hereafter LTJt) - 
with T-Type > 4, including Sc, Scd, Sd, Sdm, Sm, and Im galaxies. 

We stress that throughout this work, we do not re-classify 
local galaxies and every time we specify the morphological type, 
it refers to that from N&A classification. Here we want to measure 
a set of the most used morphological parameters in real conditions 
of the local sample and in the simulated conditions of non-local 
surveys in order to evaluate how the spatial resolution and data 
depth affect each parameter. ET, LT_et, and LTJt galaxies represent 
45%, 38%, and 17% of the selected local sample, respectively, 
completely consistent with the whole N&A sample. 

Irregular galaxies constitute only < 1% of all sources in the 
N&A catalogue. With such a small population, we were not able 
to provide reliable statistical analyses in comparison to the other 
two broader groups, Ell/SO and spirals. This is the reason why we 
do not study irregulars separately, and divide instead all late-type 
galaxies into two subgroups: earlier and later. In our LTJt defined 
group, galaxies classified as Irr by N&A form only 2% of the 
whole population, being therefore insignificant for affecting the 
results represented in this paper. 


2.2 Non-local samples 

In this Section we describe the properties of non-local surveys used 
for studying the impact from spatial resolution and data depth on 
the morphological classification of galaxies. Since our work is sta¬ 
tistical, to select the surveys we considered the following criteria: 

- a large number of detected galaxies, either a deep survey 


3000 tested local galaxies 





T-Type (N&Aclassificalion) 

Figure 1. Redshift (top), g band magnitude (middle), and morphological 
type (bottom) distributions of the local sample used in this work. 

(A colour version of this figure is available in the on-line journal) 

covering a smaller area, and/or a shallower survey covering a large 
area, 

- surveys with accurate photometry as well as accurate photo¬ 
metric redshift measurements available, and 

- surveys with different properties in terms of spatial resolu¬ 
tion and depth. 

Einally, as a result of these criteria, we selected three non-local 
surveys for our analysis: deep ground-based surveys ALHAMBRA 
and SXDS (ALHAMBRA being shallower than SXDS, but with a 
large covered area and accurate photometric and photometric red¬ 
shift measurements, and SXDS as an example of the deepest avail¬ 
able ground-based data), and COSMOS as an example of deep 
space-based surveys. We had to take into account the number of 
selected surveys, since for each of them the applied computational 
procedure is time consuming. In the following sections we describe 
each of these surveys. Their basic properties are summarised in Ta- 
ble[T] We can not compare directly the results obtained between the 
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three analysed non-local surveys: first, because they have different 
depth, so the analysed magnitude cuts are different (see below), and 
secondly they have different photometric information available, so 
the bands that we used are different as well. We tested each survey 
separately, and tried to represent the obtained results depending on 
the properties of each one. 


2.2.1 ALHAMBRA 

The Advanced Large Homogeneous Area Medium Ban d Redshift 
Astronomical survey (ALHAMBRA; iMoles et aDl2008h is a pho¬ 
tometric survey that imaged ~ 3.5 deg- of the sky, in eight different 
fields, through 20 optical and 3 near-infrared (NIR) filters. The av¬ 
erage spatial resolution of the ALHAMBRA images is ~ 1 arcsec 
(not higher than 1.6, and ranging m ainly between 0.8 and 1.2 arcsec 
in the F613W band; see Table 1 in IPovic et alj|2013l , for more de¬ 
tails). More than half million sources were detected down to the 
magnitude limit r = 25.0. 

The Bayesian Photometric Redshift code (BPZ2. Q) was 
used to measure t he photometric redshifts l lBemtez et al.l l200(il ; 
iMolino et alj|2014h . obtaining the accuracy for galaxies brighter 
than 22.5 and with magnitudes 22.5 - 24.5 in the constructed 
F814W band of ~ fiz /(1 -H z) = 0.011 and ~ fiz/(l -H z) = 0.014, re¬ 
spectively. All images and catalogues are available through the AL¬ 
HAMBRA webpag^l. 


ALHAMBRA sample selection 

In th is work we used the same sample as described in IPovic et akl 
( l2013h . whereby we will provide here only a brief summary. We 
selected only sources with the BPZ quality parameter ODDS > 0.2, 
expecting no more than 3% of photometric redshift outliers 
jMolino e t al.ll201 4l). For details about the galaxy-star separation 


Povic et alT i 2013h . And finally, we selected only objects with 


magnitudes <23.0 in the F613W band, taking into account that 
above this magnitude limit S/N ratio decreases significantly, lead¬ 
ing to less accurate photometric redshift estimations and geomet¬ 
rical galaxy-star classifications. The F613W photometric band was 
selected due to it s higher signal-to-noise (S/N) ratio in comparison 
with other filters jAparicio- Villegas et al.l[20l^ . The final selected 
sample has more than 43,000 galaxies. 


2.2.2 SXDS 

The Subaru/XMM-Vewto« Deep Survey (SXDS; TSekiguchi et al.l 
l2004h is a large, multiwavelength, and one of the deepest op¬ 
tical ground based surveys, cobvering five continuous, rectangle 
subfields, and a total area of 1.22 deg^. All subfields were ob¬ 
served in five broadband filters BVRJ'z', detecting in each more 
than 800,000 sources down to the limiting magnitudes B = 28.4, 
V = 27.8, Rc = 27.7, i' = 27.7, and z' = 26.6 at 3cr and a typical see- 
in g of 0.8. All ima ges and photometric catalogues are described 
in lFurusaw^ ( l2008h . and can be downloaded from the SXDS web- 
pag^E 

The photometric redshift distribution of each sample of galax¬ 
ies was determined using mu lti-band photometry of objects in the 
K-band-selected catalogue of [Hartley et al.l ( l2013li . having a depth 


of K = 24.3. In addition to the JHK photometry, photometric mea¬ 
surements were obtained for the uBV RJ'z! and Spitzer[fKAC chan¬ 
nels 1 and 2 using the method of ISimnson et al. I (l20lj). Pho¬ 
tometric redshift s and object classific ations were derived using 
the code EAzY jBrammer et al. 200^ using a development of 
the method described in ISimpson et all ( 1201 3h and Caldwell et 
al. (in prep). Measured photometric redshifts have accuracy of 
(5z/(l -hz) = 0.03L 


SXDS sample selection 

The sample we selected in the SXDS field is based on the available 
photometric redshifts. Therefore, we used only those galaxies that 
overlap with the K-band survey, instead of the entire SXDS field. 
The survey covers 0.77 deg^, where the overlap region with the 
SXDS field is ~ 0.63 deg^, covering the entire central subfield and 
fractions of the side ones. We discarded all galaxies with poor pho¬ 
tometric redshift fits, having < 100. To separate between point¬ 
like an d extended sources, w e applied the classification described 
above iSimpson et alJiToTsL and Caldwell et al. in prep). More¬ 
over, we selected all galaxies with i' < 24.5, which is the band that 
we will use in the following analysis, and the limit at which our 
sample is complete, taking into account the available redshifts and 
photometry. The final selected sample consists on more than 68,000 
galaxies. Since the photometric redshift selection is based on the 
K-band, we might be missing the bluest sources detected in the 
SXDS. We did not perform a detailed analysis with respect to this, 
but we estimated by cross-matching i- and K-band catalogues an 
upper limit of < 9% of these sources, which should not affect sig¬ 
nificantly our results. 


2.2.3 COSMOS 


The Cosmic Evolution Survey (COSMOS; IScoville et al.ir2007h is 
a deep, 2 deg^ multiwavelength survey. In this work we use the ob¬ 
servations from the HST ACS survey, where more than 1.2 x 10* 
sources were detected down to a limiting magnitude 26.5 in the fil- 
ter F8141V. The ACS d ata reduction and images are described in 
I Koekemoer^taL (2007), while the photometric catalogue is pre- 
sented in lLeauthaud et alj ( l2007h . 

The photome tric redshifts used in this work are presented in 
lllbert et ^ ( l2009l) . and were measured using the Le Phare code and 
photometric information from 30 broad, intermediate, and narrow- 
band filters from UV, optical, NIR and MIR bands. The obtained 
measurements show accuracy of ~(Sz/(I -(-z) = 0.007 for galaxies 
brighter than P - 22.5. At magnitudes fainter than P < 24.0 the 
accuracy is ~5z/(l-l-z) = 0.012, while for the very faint sample 
(P < 25.0) the accuracy drops significantly with ~ 20% of outliers. 
All images and catalogues used in this work are available from the 
COSMOS webpagtH 


COSMOS sample selection 

The COSMOS 2deg^ field is covered with 81 HST/ACS im¬ 
ages jKoekemoer et alhoOTl) . Taking into account the number and 
size (20,480 pixels square) of the images, and the computing time 
that we would need to process all the data with the galSVM code 
(SecO, we are finally using an area of 0.5deg“. Considering the 
depth of the COSMOS data, for the study presented in this pa¬ 
per, the selected area provides us with the statistical information 
completely comparable with the other two surveys. Moreover, we 


^ http://alhambrasurvey.com 

^ http://www.naoj.org/Science/SubaruProject/SXDS/index.html 


http://cosmos.astro.caltech.edu 



























































Morphological classification of galaxies: impact from survey depth and resolution 5 


Table 1. Summary of the ground and space-based high-redshift data selected for studying the observational bias of morphological parameters. 


Survey 

tot. A/used A [deg-] 

Sp. Res. [arcsec] 

Lim mag/Sel. sam. (band) 

Num of gal 

mag sim 

z sim 

ALHAMBRA 

4.0/~3.0 

1.1 

25.0/23.0 (F613W) 

>43,000 

< 20 . 0 , < 21 . 5 , <23.0 

- 0 . 45 ,- 0 . 62 ,- 1.0 

SXDS 

1.22/-0.65 

0.8 

27.7/24.5 (/') 

> 68,000 

<21.0, <23.0, <24.5 

-0.65,-1.2,-2.2 

COSMOS 

2.0/-0.5 

0.09 

26.5/24.0 (F8141T) 

> 107,000 

<21.0, <23.0, <24.0 

-0.67,-1.1,-1.6 


* Column description: Survey - selected high-redshift surveys; tot. A/used A - the total covered area of each survey, and the area that covers the data used 
in this work; Sp. Res. - the average survey spatial resolution; Lim. mag / Sel. sam. (band) - limiting survey magnitude, and limiting magnitude of the 
selected sample in the correspondent band; Num. of gal - the final number of selected galaxies used in this work to obtain the magnitude and redshift 
distributions in each survey, and move later the local galaxies to these conditions (see Sec.[3]and |4); mag_sim and z_sim - magnitude cuts and the 
corresponding redshifts (for 95% of the sample) used to simulate the local galaxies; 


checked that the magnitude and redshift distributions of galaxies in 
the selected field correspond to the total ones. Our main criteria for 
the COSMOS sample selection is again based on the accuracy of 
the photometric-redshifts. We selected only sources with t < 24.0, 
since the photo-z accuracy decreases significantly at fainter mag¬ 
nitudes, as described above. For the ga laxy-star separat i on we 
use d the classification s presented in both iLeauthaud et~^ 1 1200^ 
and lllbert et al.l ( l2009h . Finally, the selected sample has more than 
107,000 galaxies. 


3 METHODOLOGY 

Non-parametric methods of galaxy classification are usually ap¬ 
plied either when dealing with high redshift and/or low resolution 
data, where the galaxy decomposition and profile fitting becomes 
impossible, or in large surveys, where due to the high number of 
detected sources it is necessary to use more automated and statisti¬ 
cally based approaches. We test how much the morphological pa¬ 
rameters described in Sec. [3d] (basically related with the distribu¬ 
tion of light within the galaxy and its shape) are sensitive to depth 
and spatial resolution, and how this introduced observational bias 
might affect the final galaxy classification. To do this, we used a 
sample of local galaxies with available detailed visual classifica¬ 
tion (see Sec. ED. and measured their morphological parameters 
in two cases: 

1) at their real redshifts (magnitudes), i.e. at z~0 (see 
sec.lra. and 

2) moving them to higher redshifts (therefore to fainter 
magnitudes) and lower resolution, to simulate the conditions of 
galaxies on deep, ground- and space-based non-local surveys (see 
Sec.ED. 

In both cases all parameters are measured in a completely con¬ 
sistent way, using their definitions as described in the following 
section. Finally, we compared the results obtained in cases 1) and 
2) to quantify the impact from data depth and spatial resolution on 
each morphological parameter, for each of the three analysed sur¬ 
veys (see Sec. l4.2l and 14.3b . 


3.1 Description of the tested morphological parameters 

The six morphological parameters described here are by 
large, the most commonly used to distinguish between ET and 
LT galaxies, and to classify perturbed galaxies and interact¬ 
ing/merging systems. In the following definitions, when necessary, 
the ga laxy centre is determ ined by minimizing the asymmetry 
index dAbraham et al.l [79961 see below), while the total flux is 
defined as the one contained within 1.5 times the Petrosian radius 
(ip, measured by SExtractor; see Huertas-Company et al. 2008). 


- Abrah a m c oncentration index (hereafter CABR; 
[Abraham et alj 1 19961) - measured as the ratio between the 
flux at 30% of the Petrosian radius (F30) and the total flux (E,;,,) : 


CABR = F3 o/F,„,. 


( 1 ) 


It gives values from 0 to 1, where higher values correspond to 

higher central light fractions. _ 

- Gini co efficient (hereafter GINI; [Abraham et al.l [20031 ; 
[Lotz et si .|2004[) - presented as the cumulative distribution function 
of galaxy pixel i values: 


GINI = 


1 

\X\n(n - 1) 


Yp-i -n- DKl, 


( 2 ) 


where n is the total number of pixels in a galaxy, X, the pixel flux 
value, and X| the mean over all the pixel flux values, using as an 
aperture 1.5 r,,. Usually, it correlates with CABR, having also val¬ 
ues from 0 to 1, where higher ones correspond to galaxies with 
higher central concentrations. However, unlike CABR, it can dis¬ 
tinguish between galaxies with shallow light profiles and those with 
the light concentrated in a few pixels, but outside the galaxy centre. 

- Conselice-Ber shady concentration index (hereafter CCON; 
[Bershadv et al.l200(ih - measured as the logarithm of the ratio of the 
circular radii containing 80% and 20% of the total flux: 


CCON = 5/og(r8o/r2o). (3) 

In general, lower CCON values correspond to lower fractions of 

light in the central region. _ 

- M20 moment of light (hereafter M20; [Lotz et aD[2004h - 
measured as the flux (f,) in each pixel (x,, y,) multiplied by the 
squared distance to the centre (x^., y^) of the galaxy, summed over 
the 20% brightest pixels, and normalized by the total second-order 
moment (M,„,): 

M20 = logi ), while Vf <0.2 /„„ (4) 

^tOt 

where 

n 


It shows the light distribution of central sources, off-centre star 

clusters, bars, spiral arms, etc. _ 

-Asymmetry (hereafter ASYM: [Abraham et al.[[l99^ - mea¬ 
sured by subtracting the galaxy image rotated by 180° from the 
original image, taking also into account the background: 


V(|/(i,/)-/i 80 (i,y)|/ 2 ) 

ASYM = — -=- 


Y(\ S(h7)-Si8o(h;)l/2) 


,( 6 ) 
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where is the flux in the (i,j) pixel position on the original 
image, Ii8o(i, j) is the flux in (i, j) on the image rotated by 180°, and 
and Bi8o(i, y) are fluxes of the background in the original 
and rotated image, respectively. Defined in this way, this parameter 
quantifies the degree to which the light of a galaxy is rotationally 
symmetric. It gives values from 0 to 1, where the higher ones 
correspond to more asymmetric galactic shapes. 

- Smoothness or dumpiness (hereafter SMOOTH; 
IConselice et al.l l2000h - quantifies the degree of small-scale 
structure. To measure SMOOTH, the original galaxy image is 
smoothed out with a boxcar of a given width and then subtracted 
from the original image: 


SMOOTH = ^ -=- 


^(1 B{i,j)-Bs(i,j)\l2) 


where present again the flux (background) in the 

pixel position (i,J) on the original image, while 
Bs(i,j) ars the flux and the background, respectively, on the 
image smoothed by a boxcar of width 0.25rp. The residual image 
provides with information about possible regions of dumpiness, 
quantifying the level of small-scale structures. 


3.2 Methodology applied on the local sample 

We measured the reference morphological parameters of the lo¬ 
cal sample, at their real redshifts z~0, in three different SDSS 
bands: g, r, and i. This is needed in order to compare in a con¬ 
sistent way the morphologies of local and simulated galaxies in 
Sec. 14.21 and 14.31 since the corresponding k-correction will im¬ 
ply a different filter for such a comparison. We used the function 
implemented in the galSV]v(f| public code of galaxy classification 
dHuertas-Companv et alJl^OSl , 1200^ for measuring the parame¬ 
ters. The input information that the code needs, related with the 
source position, size, ellipticity parameters, and background, was 
obtained by running SExtractor in all three bands. 


3.3 Methodology applied at higher redshifts 

To simulate the conditions of non-l ocal surveys (case 2 in Sec.[^ , 
we u sed again the galSVM code jHuertas-Companv et aLll^OSl 
l2009li . performing the following three steps for each of the three 
analysed surveys: 

- First, we constructed the correspondent magnitude and 
redshift distributions of non-local galaxies (see Sec. 14.21 and 
Fig. H). We randomly redshifted and scaled in luminosity the 
selected sample of local galaxies to match the distributions of 
non-local ones (see Tab. [T] fot the used photometric band in each 
survey). In this way, when scaling the galaxy in flux, the surface 
brightness dimming was directly taken into account. Moreover, we 
re-sampled the local galaxies with the corresponding pixel-scale 
for each selected non-local survey, and convolved with its PSF to 
match the spatial resolution. 

- Second, we dropped the simulated galaxies, obtained in 
the first step, into the real background of the high-redshift survey 
images. We made sure that the galaxies are dropped into empty 
regions of sky, minimising the chance of superposition with 
high-redshift foreground/background sources. To detect the empty 


sky regions, the code makes use of corresponding SExtractor 
segmentation images. Moreover, although the probability of 
overlapping is related with the survey properties, we’ve been strict, 
dropping only 10 galaxies per image. This step is then repeated 
until dropping all simulated galaxies, generating every time new 
images. With all this, we make the probability of overlapping and 
merger confusion practi c ally n egligible (for more information see 
IHuertas-Companv et al.l ( 1200 sh ). Finally, by dropping the local 
galaxy in a real background of the high-redshift sources, we expect 
to reproduce the noise from the real images. 

- Third, we measured the morphological parameters described 
in Sec. 13.11 We took care of the k-correction effect introduced 
by cosmological redshift, and depending on the band selected in 
^.^^each survey, and also depending on the redshift to which the local 
galaxy was shifted, we measured the morphological parameters 
using the corresponding SDSS rest-frame band image. In the 
ALHAMBRA survey, since we are using the F613W AUTO mag¬ 
nitudes, the SDSS r band images were used for galaxies simulated 
to z_sim< 0.13, g band for the redshift range 0.13 - 0.5, and u band 
for higher redshifts. In the SXDS survey, using the i' magnitudes, 
the r band was used for redshifts z_sim<0.4, the g band for 
0.4 < z_sim < 0.86, and the u band for z_sim>0.86. Finally, in 
the COSMOS survey, using the F814W AUTO magnitudes, the i 
band was used for the lowest redshifts z_sim <0.18, the r band for 
0.18 < z_sim< 0.5, the g band for 0.5 < z_sim < 1.0, and the u band 
for z_sim>1.0. Since the filter efficiency for the SDSS u ban d 
is much lower than that for g, r, and i bands dOunn et all 1 19^ . 
and taking into account that we are analysing three broad mor¬ 
phological groups (ET, LT_et, and LTJt) as described in Sec. 12.11 
(instead of studying a finer morphology), we used the g band 
images instead of u when necessary. Figure shows an example 
of u, g, r, and i images of four galaxies from our local sample 
with different morphologies; note the poor information available 
in the u image in comparison to the other bands. Finally, since 
the morphological classification directly depends on the source 
brightness, in each survey we analysed the observational bias at the 
three different magnitude cuts: magUcut < 20.0, mag2_cut < 21.5, 
and mag3_cut < 23.0 in the F613W band in ALHAMBRA; < 21.0, 
<23.0, and <24.5, in the i' band in SXDS; and <21.0, <23.0, 
and <24.0, in the F814W band in COSMOS. In all surveys we 
chose the first magnitude cut such that we have a sufficient number 
of sources to perform our analysis. The cut at which our selected 
photometric/photometric redshift sample is complete, is the last 
one. Finally, we chose a third cut intermediate to these two. 

Figures |3]|4l and show examples of local galaxies, with dif¬ 
ferent morphologies, after being scaled to the conditions of the AL¬ 
HAMBRA, SXDS, and COSMOS surveys, respectively. In each 
survey, we show the galaxies being redshifted to the corresponding 
magnitude cuts, as explained above, providing for each cut the sim¬ 
ulated values of magnitude and redshift. These images represent the 
output of the first step explained in this section, that we use in the 
second step to drop them in the real background of the non-local 
survey images, and to measure their morphological parameters in 
the third step. The colour images wer e gathered from the EFIofl 
project, by combining the gri images dBaillard et al.ll201 ih . while 
the original images were downloaded from the SDSS DR4 database 
and correspond to the images used by N&A in their visual classifi- 


^ http://gepicom04.obspm.fr/galSVM/Home.html 


® Extraction de Formes Idealises de Galaxies en fmagerie; 
http://www.astromatic.net/projects/efigi 
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Figure 2. An example of Ell, Sa, Sc, and Sdm galaxies (from top to bottom) 
seen in four SDSS bands, u, g, r, and i (from left to right). 


cation. We can observe in each survey how the galaxy information 
changes when going from brighter to fainter magnitudes (in general 
from lower to higher redshifts), but even more, how the role that the 
spatial resolution plays when classifying galaxies, where in COS¬ 
MOS survey the galaxy information can be conserved up to much 
higher redshifts in comparison with the ground-based surveys. We 
studied this in detail in Sec.|4]and|^ 


Noise ejfect 

With respect to the noise, since we are using very bright local 
galaxies, galSVM assumes that the noise from the galaxies them¬ 
selves is negligible. Therefore, the noise coming from the galaxy it¬ 
self once it has been dropped into the non-local image is not treated 
by the code. We tested how strong this effect could be on the mea¬ 
sured morphological parameters. We selected 30 galaxies, 10 from 
each morphological group defined in Sec. irn and we measured the 
morphological parameters in the case of ALHAMBRA and COS¬ 
MOS in two cases: 1) without the above noise added, like in our 
analysis, and 2) adding the random Poisson noise, typical of each 
survey, to the scaled local images (before dropping them into the 
background of non-local surveys). We then compared the results 
obtained in these two cases for the whole samples, and in relation 
with the morphological types. We found insignificant differences of 
the above noise on the mean surface brightness (< 1%), and on the 
parameters CABR, CCON, GINI, and M20 (< 6% and < 5% in AL¬ 
HAMBRA and COSMOS, respectively, at the faintest magnitude 
cuts), but more significant differences in the case of ASYM and 
SMOOTH. In the case of ALHAMBRA, the differences for these 
two parameters go up to ~ 40% and 50%, respectively. In the case 
of COSMOS, differences for ASYM are 13%, 19%, and 57% at 
F814W < 21.0, < 23.0, and < 24.0, respectively. For SMOOTH the 
corresponding values are 15%, 35%, and 50%. From this result, 
we suggest that for the parameters ASYM and SMOOTH the noise 
effect should be taken into account when using them in morpho¬ 
logical classification of galaxies, and especially when dealing with 
ground-based data. Therefore, in our following analysis we discuss 
all results only for the four parameters CABR, GINI, CCON, and 
M20, that are more stable to noise and where our galS VM measure¬ 
ments are reliable. The results based on ASYM and SMOOTH are 
provided only in the case of COSMOS, for the first two magnitude 
cuts, where the noise effect is lower. 


4 RESULTS 

4.1 Morphological parameters of the local sample at z ~ 0 

For each galaxy, we measured the morphological parameters de¬ 
fined in Sec. 13.11 in the g, r, and i bands, as was described in 
Sec. 13.21 To test how sensitive are morphological parameters to the 
used wavelengths and if somehow it could affect our further anal¬ 
ysis, for each parameter, we compared its measurements obtained 
in the three photometric bands, finding, in general, a good linear 
correlation in all cases. Figure |6] shows the comparison of all six 
parameters between g and r (panel a), and g and i bands (panel 
b). In both cases, the concentration parameters, CABR, CCON, 
and GINI, show the best linear correlations, having similar Pear¬ 
son Correlation Coefficient (PCCfl ~ 0.9. For M20, ASYM, and 
SMOOTH, the linear correlation is also conserved in both cases, 
however they are more sensitive to the selected band in compar¬ 
ison with the concentration indexes, and show lower PCC values 
(between ~ 0.7 and ~ 0.85). They also show a slightly higher dis¬ 
persion in panel b, in comparison with panel a, with a difference in 
the PCC of ~ 0.05. M20 is the most sensitive to the selected band, 
with PCC = 0.733 between the g and r bands, and PCC = 0.688 be¬ 
tween the g and i measurements. The number of catastrophic out¬ 
liers, placed in the ’wings’ above and below the linear correlations 
is ~ 10% in both panels. We checked the properties of these sources 
(magnitude, size, and redshift), without finding any significant dif¬ 
ference in comparison with the whole sample. We only found a 
small difference in relation with the morphological type, with 6% 
more LTs in the wings. To have an idea about how these outliers 
could affect the results in our following analysis and if there are 
some important differences that we should study in more detail, we 
performed a test by excluding them from the whole sample and 
checking again the corresponding distributions of morphological 
parameters (see Sec. 14.21 and Figs. [S] and llOt . without finding 
any significant differences. However, this was done only to test 
the properties of outliers, but the all analysis and results presented 
in this paper were obtained using the whole sample, outliers in¬ 
cluded. For the comparison of mo rphological parameters betwee n 
near-infrared and optical bands see iHuertas-Companv et al.l(l2014l) . 


4.2 Morphological parameters of the simulated sample 

To measure the morphological parameters of the simulated sample, 
we followed the methodology described in Sec. 13.31 We first 
moved randomly the local sample to map the magnitude and 
redshift distributions of galaxies in each high-redshift survey, at 
each selected magnitude cut. Figure |3 shows the corresponding 
magnitude and redshift distributions of the simulated galaxies in 
case of ALHAMBRA (panel a), SXDS (panel b), and COSMOS 
(panel c). Since the magnitude cut introduces also a redshift selec¬ 
tion, we are dealing at the same time with different redshift ranges 
in the three surveys (as summarized in Table [T). After simulating 
the local galaxies and dropping them into the real background 
of high-redshift surveys, we run again SExtractor to measure 
the same parameters as we did at z ~ 0. Finally, we measured 
the morphological parameters of simulated galaxies, again in the 
same way as we did at z ~ 0. As mentioned above, in each survey, 
depending on the used photometric band and the redshift to which 


^ PCC is a measure of the linear correlation between two variables giv¬ 
ing a value -1 <PCC < +1, where -l-l shows total positive c orrelation, 0 no 
correlation, and -1 total negative correlation jPearso 
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Figure 3. Example of simulated images of local galaxies (with different morphologies) after being scaled to map the conditions of the ALHAMBRA survey. 
For each galaxy (in each row) we present the following (from left to right): colour image, used rest-frame image(s), and simulated high-redshift images at 
three magnitude cuts (as written on the top of each redshifted image). The morphological type of each galaxy is noted on the colour image. The real redshift 
and magnitude of the galaxy are noted on the corresponding rest-frame band image(s), while the simulated high-redshift and magnitude are noted on the 
scaled images of each magnitude cut. 


(A colour version of this figure is available in the on-line journal) 


the galaxy was moved, we used the corresponding rest-frame 
SDSS images to carry out all the measurements. 

In each non-local survey, we had to repeat all the procedure 
a number of times, for each separate image, in order to cover the 
analysed survey areas (see Table [T). In case of ALHAMBRA, we 
were dealing with 48 sub-fields (4 CCD images per pointing, wit h 
12 pointing in total (1 or 2 per each field); see lPovic et alJuOl^ . 
In SXDS we repeated the procedure in each of five SXDS fields 
(see Sec. I2.2.2t . while fo cover fhe selected area of COSMOS 
we run fhe code on 17 images. Moreover, in all surveys we 
repeated fhe measuremenfs three times, since having selected three 
magnitude cuts (see Table [T). Since each time we run galSVM 
on one image, 3000 simulated galaxies were distributed randomly 
to the corresponding redshift and magnitude distributions of 
the non-local samples, for each galaxy we obtained a number 
of measurements that is equal to the number of used images. 
Therefore, the final catalogues contain 144000, 15000, and 51000 


simulated galaxies in the case of ALHAMBRA, SXDS, and 
COSMOS, respectively. Each catalogue contains the columns 
with morphological parameters measured at higher-redshifts in 
the three magnitude cuts, besides the corresponding reference 
measurements (obtained at z ~ 0). Once again, when assigning 
the reference value to the simulated one, we took care of the 
k-correction, assigning the value measured in the corresponding 
rest-frame band (that will depend on the photometric band we 
used for analysis in each non-local survey, and on the simulated 
redshift). The reader should keep this in mind for understanding 
easier the distributions of comparison sample of local galaxies, 
represented in figures from[TT]to[T7]in Sec. 14.31 

The obtained morphological parameters of simulated galaxies 
are represented with the filled-line histograms in Fig. and 
HOI that reproduce the conditions of the ALHAMBRA, SXDS, 
and COSMOS surveys, respectively. In each figure are shown the 
distributions of CABR, GINI, and CCON concentration indexes, 
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Figure 4. Same as Figure|3] but in the SXDS survey.(A colour version of this figure is available in the on-line journal) 


M20, asymmetry, and smoothness (from left to right panels) for 
three morphological groups: ET (top plots), LT_et (middle plots), 
and LTJt (bottom plots), defined in Sec. 12.11 For each morpholog¬ 
ical group, we represent the distribution of parameters at the three 
analysed magnitude cuts. In these figures we also represent the 
distributions of the reference morphological parameters of local 
galaxies, measured in Sec. 14. II with the filled-colour histograms. 

For all morphological parameters we can observe important 


differences between the reference (local) values and those obtained 
after moving the local galaxies to the conditions of non-local sur¬ 
veys, indicating the influence that the observational bias from spa¬ 
tial resolution and depth can have on them. To quantify this bias, 
we measured the difference between the median values of the local 
and simulated distributions, normalised with respect to the median 
value of the local sample: 

< local > - < simulated > / < local >. 
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Figure 5. Same as Figure[3] but in the COSMOS survey.(A colour version of this figure is available in the on-line journal) 


The obtained values are expressed in % in each plot, together 
with the corresponding range of interquartile range (IQR; between 
parentheses), measured as the difference between the third and first 
quartiles. In the following, we describe the obtained differences, 
considering as significant those within > 2 IQR. 

CABR 

In ALHAMBRA conditions, CABR showed significant differ¬ 
ences in the case of the faintest ETs and brightest LTs (~ 30% 
within 2 IRQ), with galaxies appearing less and more concentrated, 
respectively. In the conditions of SXDS we did not find any dif¬ 


ferences (all within 1 IQR). In COSMOS, we obtained differences 
of ~ 20 - 36% (within 2 IQR) at the brightest magnitude cuts for all 
three morphological types, with simulated galaxies appearing more 
concentrated. In ALHAMBRA, 2-3% of sources showed invalid 
measurements (with values outside the typical range of this param¬ 
eter, between 0 and 1), independently of the morphological type; in 
SXDS, this went between 3 - 5% for ET and LT_et, and 7 -10% for 
LTJt galaxies; while in COSMOS, 2-3% for ET and LT_et, and 
8 -10% for LTJt galaxies. 
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Figure 6. Comparison between the six morphological parameters for the 
galaxies in the local sample measured in the g and r bands (panel a), and g 
and i bands (panel b). In each panel we compared, from top to bottom, and 
from left to right: Abraham concentration index, asymmetry parameter, 
Gini coefficient, smoothness, Conselice-Bershady concentration index, and 
M20 moment of light. In all plots the red solid line presents the perfect 
linear congelation between the two measurements. 

(A colour version of this figure is available in the on-line journal) 


GINI 

In ALHAMBRA, this parameter showed differences of 10-20% 
(2IQR) at the first magnitude cut for all types, and at the second 
magnitude cut for LT_lt galaxies. In SXDS, the difference for LT_et 
galaxies was observed at all three magnitude cuts (15 -18% within 
2 IQR). The same was obtained for LTJt sources but with slightly 
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Figure 7. Magnitude (left plots) and redshift (right plots) distributions of 
simulated galaxies for the conditions in ALHAMBRA (panel a), SXDS 
(panel b), and COSMOS (panel c) surveys. In each survey, they are rep¬ 
resented for three magnitude cuts, from top to bottom: for F613W <20.0, 
F613W<21.5, and F613W<23.0 in the ALHAMBRA; /<21.0, i<23.0, 
and i < 24.5 in the SXDS; and F814 < 21.0 (top), F814 < 23.0 (middle), and 
F814 < 24.0 in the COSMOS. 
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higher differences. In the case of COSMOS, the only differences 
were observed at the first magnitude cut: ~ 20% (2 -3IQR) for all 
three types. In all cases where the differences were detected, we ob¬ 
tained negative values, showing that the simulated galaxies would 
appear more concentrated than they really are. With respect to the 
measured values of simulated galaxies that are outside of the valid 
ranges, the populations are in agreement with those for the CABR 
parameter. 

CCON 

A different behavior of CCON was detected for the conditions 
of the analysed ground- and space-based surveys. In the case of 
COSMOS, we did not detect any differences, independently of the 
magnitude cut and morphological type. This index showed to be es¬ 
pecially sensitive to the spatial resolution in comparison to CABR 
and GINI. In the case of ALHAMBRA and SXDS, the highest dif¬ 
ferences are detected for ET galaxies at the highest magnitude cuts, 
~ 30% at 3 IQR and 4 IQR. In all cases where differences were ob¬ 
served in the ground-based surveys, lower concentrations appeared 
for simulated galaxies in comparison with the reference (local) val¬ 
ues, independently on the magnitude cut and morphological type. 
Again, the population of the obtained invalid measurements is con¬ 
sistent with those presented for CABR. 

M20 

As described in Sec. 13.11 this parameter is sensitive to galaxy 
structures as the central light from the bulge, bars, spiral arms, 
strong star formation regions, etc. Therefore, for local ET galax¬ 
ies it mainly maps the central light, having a much narrower distri¬ 
bution than in the case of LT_et, or even more in the case of LTJt 
galaxies, where the distribution is broader (filled-colour histograms 
in Figs.l^O and llOb . However, when moved to higher-redshift sur¬ 
veys the distribution of this parameter changed significantly (both, 
shape and range; filled-line histograms), and showed in all cases the 
information coming mainly from the centre of the galaxy, and inde¬ 
pendently on the morphological type. In ALHAMBRA, the highest 
differences are of order > 30% within 2 -3 IQR, for all three types. 
In SXDS conditions, for ET galaxies M20 suffered similar changes 
as in ALHAMBRA (at higher magnitude cuts), except for LT_LT 
sources where no difference were found respect to local sample. In 
COSMOS conditions, M20 shows higher concentrations for sim¬ 
ulated galaxies in comparison with the local ones, affecting only 
LTJt galaxies, with differences <-31% (within 2 IQR). The popu¬ 
lation of sources with invalid measurements of M20 is consistent 
with that for CABR. 

ASYM 

This parameter showed the highest differences in the case of AL¬ 
HAMBRA (even up to ~ 80% within 3 IQR at the faintest magni¬ 
tude cut F6I3W < 23.0), affecting LT galaxies in the sense of ob¬ 
taining more symmetric galaxies in comparison with their reference 
(z ~ 0) values. In COSMOS, all three morphological types suffered 
similar shifts (~ 30% within 2 IQR), being the strongest one at the 
lowest magnitude cut F814W <21.0, but showing positive values in 
comparison to ALHAMBRA, and so producing more asymmetric 
galaxies in simulated conditions. The population with invalid mea¬ 
surements was higher (<20%) than in the case of CABR, GINI, 
and CCON. All these measurements were obtained without tak¬ 
ing into account the noise introduced by the galaxy itself once be¬ 
ing dropped into the background of non-local surveys. However, 
as showed in Sec. 13.31 in the case of ASYM and SMOOTH this 
noise could have a significant effect on these parameters, and the 


differences measured there should be added to those measured in 
this section. Therefore, we found these parameters unreliable for 
the use at higher redshifts in the case of ground-based surveys. In 
the space-based surveys similar to COSMOS, it becomes unreliable 
for the use in the case of galaxies fainter than F814W = 24.0. 

SMOOTH 

Without taking into account the noise studied in Sec. 13.31 the dif¬ 
ferences in SMOOTH in the case of ALHAMBRA were detected 
only for LT galaxies, being lower than 36% (within 2 IQR); in the 
case of SXDS, no differences were observed between the simulated 
and local samples; nevertheless, in COSMOS conditions it showed 
40 - 50% differences (within 2-3 IQR) at the brightest magnitude 
cuts F8I4W<2L0, resulting in galaxies to appear more clumpy. 
Again, to these differences those related to the inclusion of noise 
(see Sec. should be added. Moreover, this parameter showed 
to be significantly more affected with both spatial resolution and 
depth, when considering the number of invalid measurements of 
the simulated sample (up to 50% and 20% at the faintest magnitude 
cuts in the case of ALHAMBRA and COSMOS, respectively). 
Taking all this into account, we found this parameter especially 
sensitive to all, noise, spatial resolution, and data depth, and 
therefore unreliable for the use in the morphological classification 
of galaxies in the case of ground-based surveys, and for faintest 
sources in space-based surveys. 


4.3 Behaviour of morphological diagnostic diagrams 

In this section we analyse how the observational bias, discussed in 
the previous section, affects some of the commonly used diagnostic 
diagrams, that allow to separate between the ET and LT galaxies, 
to select disturbed galaxies, and/or interacting (merging) systems. 
To do this, we analysed the simulated sample in the conditions 
of the three non-local surveys at different magnitude cuts, and 
compared to the corresponding local sample. We can directly use 
the obtained diagrams to quantify how strong is the bias for each of 
these diagrams when taking into account the survey observational 
properties; basically the spatial resolution (comparing the results 
between the ground-based and space-based surveys) and survey 
depth (observing the changes at different magnitude cuts). 

As already mentioned in Sec. 13.31 we analysed diagrams 
based on CABR, GINI, CCON, and M20 in all surveys, since these 
parameters are more stable to the noise effects. Figures [TT] [T^ 
[13] [T4| [T^ and [T^show the relation between CCON and CABR, 
CCON and GINI, CCON and M20, GINI and CABR, M20 and 
CABR, and finally M20 and GINI, respectively. In each figure are 
shown the results obtained for the conditions in the three surveys: 
ALHAMBRA (top), SXDS (middle), and COSMOS (bottom). For 
each survey we plot the locus of the three analysed morphological 
groups, ET (red contours), LT_et (blue contours), and LTJt (green 
contours), obtained once local galaxies (bottom plots) were moved 
to higher-redshift conditions (top plots). Although initially we 
simulated the same sample of 3000 local galaxies to the conditions 
of the three non-local surveys, in some plots of the reference (local) 
values small differences can be observed, which is the reason why 
we represent them in each survey and for each magnitude cut. 
To understand these differences it has to be considered that, first, 
the local measurements correspond to the rest-frame band of the 
simulated sample, and hence slightly changes depending on the 
magnitude cut and survey (see Sec. 14.2t . Secondly, the simulated 
galaxies with invalid measurements of morphological parameters 
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Figure 8. Comparison between the morphological parameters of the sample of local galaxies measured at their real redshift (filled-colour histograms, repre¬ 
senting the reference galaxy properties) and in conditions that con'espond to the galaxies from the ALHAMBRA survey (filled-line histograms). Six parameters 
are compared, from left to right: Abraham concentration index, Gini coefficient, Conselice-Bershady concentration index, M20 moment of light, asymmetry, 
and smoothness. Top) three magnitude cuts of the ET galaxies, middle) LT_et galaxies, and bottom) LTJt galaxies. In all diagrams the numbers represent the 
difference between the median value of the reference (local) and simulated distributions, normalised with the median of the local sample (expressed in %), and 
the coiTesponding interquartile range. (A colour version of this figure is available in the online journal) 



















































































































































































































14 Povic et al. 


ET galaxies (SXDS conditions) 


2000 

mag_hz < 21.0 

mag 

_hz<21.0 

ma 

g_hz < 21.0 


mag_hz < 21.0 

mag_hz < 

21.0 mag_hz<21.0 

1500 

z 

1000 

500 

.± 

j 

k 


ii 

1 

lx.. 


i 

2000 

1500 

z 

1000 

500 

mag_hz < 23.0 
2%(1IQR) 

mag 

-12%(1IQF 

i 

_hz < 23.0 

) 

i 

ma 

26% (3IQR 

n ii 

g_hz < 23.0 

) 

1 

mag_hz < 23.0 

31%(3IQR) 

iL.. 

mag_hz < 
22%{1IQR) 

23.0 mag_hz < 23.0 

14%(1IQR) 

A& 

2000 

1500 

z 

1000 

500 

mag_hz < 24.5 
11%(1IQR) 

'A 

mac 

] 

_hz < 24.5 

L 

ma 

29% (4IQR 

g_hz < 24.5 

) 

1 

mag_hz < 24.5 
35% (4IQR) 

II 

mag_hz < 
17%{1IQR) 

24.5 mag_hz < 24.5 

3%(1IQR) 

k'- Ja 


0.0 0.2 0.4 0.6 0.2 0.4 0.6 0.8 1 2 3 4 -3 -2 -1 0 -3 -2 -1 0 -3 -2 -1 0 

CABR GINI CCON M20 logASYM logSMOOTH 


LT et galaxies (SXDS conditions) 


mag_hz < 21.0 

-15%{1IQR) 

A 

me 

-18%{2IQ 

i 

g_hz < 21.0 

R) 

L 

15%(2IC 

■nag_hz < 21.0 

R) 

Ifc 

mag_hz < 21.0 
19%{2IQR) 

ik 

mag_hz 

8%{1IQR) 

<21.0 

mag_hz < 21.0 

-3%(1IQR) 

mag_hz < 23.0 

-10%(1IQR) 

m 

-17% 

{2IQR) 

g_hz < 23.0 

17%(2IC 

t 

Tiag_hz < 23.0 

R) 

tk 

22% (S 

A 

mag_hz < 23.0 

IQR) 

L. 

mag_hz 

2%(1IQR) 

J: 

<23.0 

mag_hz < 23.0 

-1%(1IQR) 

mag_hz < 24.5 

-3%(1IQR) 

1 

me 

-15% 

(2IQR) 

^hz < 24.5 

_ 

20% (2IC 

J; 

Tiag_hz < 24.5 

R) 

26% (2 

A 

mag_hz < 24.5 

IQR) 

mag_hz 

-6%(1IQR) 

<24.5 

L 

mag_hz < 24.5 
-6%(1IQR) 


0.0 0.2 0.4 0.6 0.2 0.4 0.6 0.8 1 2 3 4 -3 -2 -1 0 -3 -2 -1 0 -3 -2 -1 0 

CABR GINI CCON M20 logASYM logSMOOTH 



0.0 0.2 0.4 0.6 0.2 0.4 0.6 0.8 1 2 3 4 -3 -2 -1 0 -3 -2 -1 0 -3 -2 -1 0 

CABR GINI CCON M20 logASYM logSMOOTH 


Figure 9. Same as Fig.O but moving the local sample to map the conditions of the SXDS survey. 


(A colour version of this figure is available in the online journal) 
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Figure 10. Same as Fig.[8] but moving the local sample to map the conditions of the COSMOS survey. 


(A colour version of this figure is available in the online journal) 
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(discussed in the previous section) are excluded in all plots, and 
therefore the corresponding local galaxies are also excluded from 
the bottom diagrams. ASYM and SMOOTH were only used in the 
case of COSMOS and its first two magnitude cuts. FigurefTTl shows 
some of the commonly used relations, representing the relation 
between CABR and ASYM (top plots), CCON and SMOOTH 
(middle plots), and finally ASYM and GINI (bottom plots). In all 
three plots, we represent the same comparisons as in Fig.ll lltollbl 
In all figures, when observing the reference parameters for 
local galaxies, we can clearly separate the regions typically occu¬ 
pied by ET, LT_et, and LTJt galaxies, as expected. On the other 
side, we can observe how the position and the shape of the same 
regions change once we go to fainter magnitudes (higher redshifts). 
We measured the level of contamination for the highest density 
population (50%) of each analysed morphological group (ET, 
LT_et, and LTJt) with the other two types, for the morphological 
diagrams presented in Fig. m to 1 161 We applied the Tukey’s Five 
Number Summary statistic jTukevll 1977h . and measured for each 
parameter, and each morphological type: the sample minimum, 
first quartile, median value, third quartile, and sample maximum. 
The first and third quartiles, present the median values of the lower 
half and the upper half of the sample, respectively, with respect to 
the median value of the total sample. Therefore, measuring the first 
and the third quartile of each parameter on the particular diagnostic 
diagram, we are able to determine the regions that correspond to 
the highest population (50% of sources) for ET, LT_et, and LTJt 
around their median values. From this, for each morphological 
type, we can estimate the level of contamination by the other two 
types, within their densest regions: e.g., contamination of the dens¬ 
est ET region with the densest region of the LT_et galaxies (as a 
number of LT_et galaxies in the overlapping area, normalized with 
the number of ET galaxies), and/or with the densest region of the 
LTJt galaxies, and vice-versa. Table|2]shows these measurements 
for simulated galaxies in ALHAMBRA (top table), SXDS (middle 
table), and COSMOS (bottom table), for their corresponding 
magnitude cuts. These measurements represent a quantification 
of the bias represented in the diagrams [TT] to [T^ in all the three 
surveys, and in Eig. [T^in the case of COSMOS. We can again 
observe that the contamination levels in the morphological types 
increase with the data depth in all three surveys, being however 
significantly lower in the case of COSMOS. 


5 DISCUSSION 

5.1 Individual parameters 

As seen in Sec. 14.21 the morphological parameters show complex 
behaviors when going from one survey to another, from one magni¬ 
tude cut to another, and when dealing with different morphological 
types. In principle, they suffer changes due to the impact of spatial 
resolution and survey depth, and depending on how strong these 
two effects are, they might cause: 

- only a re-distribution of the light within the galaxy, without 
changing much the total flux, resulting in a weak impact, and/or 

- the loss of the galaxy low surface brightness structures and 
the change of the total flux, hence producing a more significant im¬ 
pact. 

However, how the distributions of morphological parameters 
behave in each survey when going from brighter to fainter magni¬ 
tude cuts (from lower to higher redshifts), depends on their defi¬ 
nitions. In general, for weak impacts of spatial resolution and data 
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Figure 11. Relation between the Conselice-Bershady and Abraham con¬ 
centration indexes, in the three analysed non-local surveys: ALHAMBRA 
(top plots), SXDS (middle plots), and COSMOS (bottom plots). In each 
survey, top rows represent the morphological parameters of simulated 
sample obtained after moving the local galaxies to higher redshifts, 
considering three magnitude cuts (the first and the last column showing the 
lowest and the highest analysed magnitude cut, respectively), while bottom 
rows show the corresponding reference (local) values. The red, blue, and 
green contours represent ET, LT.et, and LTJt galaxies, respectively. 


(A colour version of this figure is available in the online journal) 
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Figure 12. Same as Fig. m\ but representing the relation between the 
Conselice-Bershady concentration index and Gini coefficient. 


Figure 13. Same as Fig. \m but representing the relation between the 
Conselice-Bershady concentration index and M20 moment of light. 


(A colour version of this figure is available in the online journal) 


(A colour version of this figure is available in the online journal) 


depth, the galaxies might appear as more concentrated, asymmetric 
and clumpy, while for stronger impacts, when the galaxy total flux 
can be affected significantly, the galaxies might appear as less con¬ 
centrated, more symmetric, and less clumpy. To trace the real im¬ 
pact from spatial resolution and depth on the parameters in Fig. [8] 
|9l and[T3 it is not only important to take into account the difference 
between the simulated and real values (and its standard deviation), 
but also the sign of that difference. In general, for the concentra¬ 
tion parameters (CABR, GINI, CCON), the distributions will be 
the most affected when the final differences have positive values 


(showing that the simulated galaxies are less concentrated in com¬ 
parison with the real ones); negative in case of M20, ASYM, and 
SMOOTH due to the use of the logarithm scale in our analysis. 

5.2 Diagnostic Diagrams: impact from spatial resolution and 
data depth 

To evaluate the separate effect of resolution and data depth, we 
compare two cases: first, the results obtained for the two ground- 
based surveys, where both of them are affected with the spatial res- 
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Figure 14. Same as Fig. ED but representing the relation between the Gini 
coefficient and Abraham concentration index. 


Figure 15. Same as Fig. \n\ but representing the relation between the M20 
moment of light and Abraham concentration index. 


(A colour version of this figure is available in the online journal) 


(A colour version of this figure is available in the online journal) 


olution, and we can study how important is the survey depth, and 
second, the results obtained for space-based data in comparison to 
the ground-based surveys, where we can see the role that play the 
spatial resolution in morphological classification. 

5.2.1 The effect of survey depth in deep ground-based surveys: 
ALHAMBRA and SXDS 

From the morphological diagrams and from the Tabled for 
the deep ground-based surveys, like ALHAMBRA and SXDS, di¬ 
agrams combined with M20 show to be less effective in separat¬ 


ing galaxies. On the other side, as seen in all figures, survey depth 
is crucial when classifying galaxies. In the case of ALHAMBRA, 
at its lowest magnitude cut (F613W < 20.0), although in all dia¬ 
grams the position/shape of locus of the three morphological types 
change significantly in comparison with the same sample at z ~ 0, 
in most of diagrams we are still able to distinguish between the 
ET and LT galaxies, or at least the regions of the highest source 
densities do not overlap. Observing for each type their contours en¬ 
closing 50% of galaxies, we found a ~ 5 - 25% contamination of 
ET with LT_et sources, and a ~ 3 - 35%, of LT_et with ET sources, 
where the contamination level depends on the used diagram, with 
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M20 M20 M20 
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M20 M20 M20 

Figure 16. Same as Fig. m but representing the relation between the M20 
moment of light and Gini coefficient. The green and red lines represent the 
limits of ILotz et al] <200 Sh to distinguish between the normal galaxies and 
mergers, and between ET and LT galaxies, respectively. In their classifica¬ 
tions, mergers occupy the regions above the green lines, ETs (LTs) below 
the green lines and above (below) the red ones. Dotted lines (both green 
and red) provided this separation for a sample of local galaxies, while solid 
lines were used to classify sources at 0.2 < z < 0.4. 

(A colour version of this figure is available in the online journal) 


GINI vs. CABR and CCON vs. CABR diagrams showing the low¬ 
est contamination (< 5%), and CCON vs. M20, and M20 vs. GINI 
the highest one. Similar conclusions are obtained for classifying 
LTs, but with higher contamination by ET sources. When includ¬ 
ing fainter galaxies (second magnitude cut F613W < 21.5), the mix¬ 
ing between the three regions increases significantly, but still for 


mag_simul < 21.0 mag_simul < 23.0 



mag_simul < 21.0 mag_simul < 23.0 


simulated galaxies 
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Figure 17. Relation between the Abraham concentration index and asym¬ 
metry (top), Conselice-Bershady concentration index and smoothness (mid¬ 
dle, top), asymmetry and smoothness (middle, bottom), and asymmetry and 
Gini coefficient (bottom) in COSMOS, in the two analysed magnitude cuts, 
where the noise effect on ASYM and SMOOTH parameters is still insignif¬ 
icant (see Sec. IS). For more explications see the caption of Fig. ED (A 
colour version of this figure is available in the online journal) 
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diagrams that do not use M20 parameter the contamination lev¬ 
els of ET and LT_et are between 20 - 30%, depending on the di¬ 
agram. In the case of the last magnitude cut (F613W<23.0) the 
mixing is much higher, ~60-100%, making each of the diagrams 
useless if used separately. On the other side, in the case of SXDS 
we observe a similar situation at the different magnitude cuts as 
in the case of ALHAMBRA, but when including fainter galax¬ 
ies at higher redshifts. The contamination levels obtained at the 
lowest magnitude cut (i' <21.0) in SXDS (taking into account the 
central wavelength of the i' band, of ~760nm, this would corre¬ 
sponds to even higher magnitudes in ALHAMBRA), are higher 
than those obtained at the lowest analysed magnitude cut in AL- 


classification of galaxies. 


5.3 Diagnostic diagrams: implications for morphological 
classification of galaxies 

For intermediate- and high-redshift galaxy samples, and large data 
sets, typical for current (and future) deep photometric surveys, 
the application of automated non-parametric methods, that use the 
parameters analysed within this work, is one of the ways to obtain 
the information about the morphology. In some of the latest works, 
the classification based on non-parametric methods was done using 


HAMBRA (F613W < 20.0), but lower than the second magnitude 
cut (F613W < 21.5). Moreover, at the highest magnitude cut in the 

a set of parameters simultaneously, therefore decreasing the impact 
of spatial resolution and depth (e.g. Huertas-Comoany et al.ll2009l 

SXDS (i'<24.5), we obtained ~ 10-80% lower contaminations, 
depending on the type, than in the case of the last magnitude cut in 

2011 

: IPovic et alJ 

20 i4 

20131; Cibineletal, 2013; Carollo et al.l 

2013 

; ICotini et al 

2013 

; iHolwerda et alj I20l4. However, for 


ALHAMBRA (F613W < 23.0). The combination of other parame¬ 
ters with M20 (Fig. [13] (Hj and El}, which seems to be the most 
affected with the resolution, again shows to be the most ineffective. 
Finally, we should notice that either the effectiveness in separating 
between ET and LT_et changes from one diagram to another and 
from one magnitude cut to another, LTJt galaxies in all cases, in 
both surveys, suffer very high levels of contamination (> 80% in 
most of cases) from other two morphological types. 


5.2.2 The effect of resolution. COSMOS as a reference 

Since the spatial resolution affects little the measurements of 
the morphological parameters in the case of COSMOS data, in 
comparison with the results obtained in ALHAMBRA and SXDS 
surveys, we can evaluate how much impact does it have when 
classifying galaxies. In all diagnostic diagrams (Fig. to llbt the 
loci of galaxies simulated to the conditions of COSMOS, irrespec¬ 
tive of their morphological type, change significantly less at all 
magnitude cuts than in case of ALH and SXDS. Moreover, the 
diagrams that include ASYM and SMOOTH parameters provide 
the additional information, and can be used up to F814W = 24.0 
magnitudes. However, the data are still affected by the survey 
depth, as can be seen in all figures when comparing the diagrams 
at different magnitude cuts. When including only the brightest 
simulated galaxies in the sample (F814W<2L0) the locus of 
the three morphological groups basically does not change in 
comparison with that of the corresponding local values, and the 
contamination levels between ET and LT galaxies are < 10%. This 
picture starts to change when including fainter galaxies, down to 
F814W<23.0, and even more at F814W<24.0. Depending on 
the diagram the contamination rate changes between 2-25% at 
F814<23.0, and 6-35% for ET and 20-70% for LT galaxies at 
F814 < 24.0. However, even at the highest analysed magnitude cut, 
which corresponds to even fainter magnitudes in ALHAMBRA 
and SXDS selected bands, we still don’t observe the complete 
overlapping of the ET and LT regions like in the case of ground- 
based surveys, and the contamination levels are still below those 
obtained for the brightest ALHAMBRA/SXDS magnitudes. 

Although the impact from survey depth on the morphological 
parameters measured within COSMOS conditions is still present, 
it is significantly lower than in the case of deep ground-based 
surveys like ALHAMBRA and SXDS. We can assume that the 
spatial resolution is the main limitation that should be taken into 
account when dealing with the photometric data and morphological 


classifications based on only one, two, or three parameters, special 
attention should be taken in relation to the magnitude and redshift 
distribution of the analysed sample, and the data quality. 

The ASYM index was used in s ome of previous wo rks to 
select mergers and disturbed galaxies ( lLotzetalJl20l03fll201ll: 
IVillforth et al.ll201^ . ASYM is especially sensitive to the survey 
noise, as showed in Sec. 13.31 and to both, the spatial resolution and 
depth, as showed in Sec. 14.21 In surveys similar to COSMOS, for 
sources fainter than F814 = 23.0, we again do not recommend the 
use of a single criterion based on ASYM when classifying them 
morphologically, since the contamination levels of 50% highest 
density regions of ET and LT rise to above 40%. In the case of 
ground-based surveys, we do not recommend the use of this single 
criterion in any case (except for maybe the brightest galaxies), 
since the contamination levels are much higher, taking the noise 
effect especially into account; we suggest instead the simultaneous 
use of multiple diagnostic diagrams. 

Over the past ten years, the M20-GINI diagram (Fig. QU 
found wide applications. It was mainly used in identifying 
galaxy mergers and peculiar (disturbed) galaxies (e.g. Lotz et al.l 


2004 pTerce et al.l l2007l; 


Lotz et; 


tur bed) galaxies (e.g. ILot zeLaU 
ZI Eooi IConselice etal.ll2008l: 
Mendez_et_alj20D]J^ sn et aD l2009l: IStott et al.ll2013l: IPettv et all 


20141 ; Hung et aH 2014 ); often it was used in the analysis where 


different types of galaxies, e.g. ULIRGs, dust obscured galaxies, 
star-forming galax ies, or quiescent galaxi e s, were compared with 
control samples ]Bussmann.eLaL 2004 Pentericci_et^ 201C 


Lanyon-Fostei^alJ |2012| ; IWang et al.l |2012| ; IChung et al.l 12013 
for selecting ET and LT galaxies jLotz et al.l 


Bdhm et al.ll2014l) 

2004l2008l) . 


From our results, in ALHAMBRA and SXDS surveys even 
for the brightest analysed sources (F613<20.0 and i'<2L0, 
respectively) the contamination levels between the ET and LT 
galaxies in M20-GINI diagram are about 30%. In COSMOS 
however, at the lowest magnitude cut F814<2L0 (z<l) the 
contamination is insignificant (< 10%, except for the youngest 
spiral/irregular galaxies), but already at the second magnitude cut 
F814 < 23.0 (z < 1.5) the contamination is > 20%. We re presen ted 
in Fig. [T^ the limits established by iLotz et all | |2004 boOSh to 
separate between normal ET and LT galaxies and major mergers, 
in the local universe (dotted lines), and at 0.2 < z < 0.4 (solid 
lines), using the All-wavelength Extended Groth strip International 
Survey (AEGIS) HST/ACS data. The solid-lines criterion (non¬ 
local sample) was later applied by iLotz et aP ( 120081 12010al f9) to 
z< 1.2 sources to select major mergers, reporting a population of 
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Table 2. Expected level of contamination for densest regions (50% of sources) of ET, LT_et, and LTJt galaxies in some of the most used morphological 
diagrams in the ALHAMBRA (top), SXDS (middle), and COSMOS (bottom) surveys. 




magl.cut < 20.0 

magl.cut < 21.5 

magS.cut < 23.0 


Diagram 

cont. ET 

cont. LT.et 

cont. LT_lt 

cont. ET 

cont. LT.et 

cont. LTJt 

cont. ET 

cont. LT.et 

cont. LTJt 



[%] 

[%] 

[%] 

[%] 

[%] 

[%] 

[%] 

[%] 

[%] 


CCON 

- (LT.et) 

4(ET) 

-(ET) 

24 (LT.et) 

31 (ET) 

2(ET) 

58 (LT.et) 

71 (ET) 

92 (ET) 

A 

vs. CABR 

3 (LTJt) 

11 (LTJt) 

42 (LT.et) 

- (LTJt) 

20 (LTJt) 

77 (LT.et) 

20 (LTJt) 

35 (LTJt) 

100 (LT.et) 

L 

CCON 

11 (LT.et) 

13 (ET) 

-(ET) 

31 (LT.et) 

38 (ET) 

6(ET) 

56 (LT.et) 

78 (ET) 

100 (ET) 

H 

vs. GINI 

- (LTJt) 

16 (LTJt) 

58 (LT.et) 

1 (LTJt) 

24 (LTJt) 

86 (LT.et) 

22 (LTJt) 

36 (LTJt) 

100 (LT.et) 

A 

CCON 

- (LT.et) 

46 (ET) 

-(ET) 

50 (LT.et) 

75 (ET) 

13 (ET) 

65 (LT.et) 

93 (ET) 

100 (ET) 

M 

vs. M20 

34 (LTJt) 

12 (LTJt) 

40 (LT.et) 

2 (LTJt) 

19 (LTJt) 

88 (LT.et) 

21 (LTJt) 

34 (LTJt) 

100 (LT.et) 

B 

GINI 

3 (LT.et) 

3 (ET) 

16 (ET) 

29 (LT.et) 

31 (ET) 

49 (ET) 

58 (LT.et) 

77 (ET) 

100 (ET) 

R 

vs. CABR 

- (LTJt) 

22 (LTJt) 

100 (LT.et) 

3 (LTJt) 

36 (LTJt) 

100 (LT.et) 

28 (LTJt) 

42 (LTJt) 

100 (LT.et) 

A 

M20 

6 (LT.et) 

8 (ET) 

-(ET) 

37 (LT.et) 

48 (ET) 

23 (ET) 

62 (LT.et) 

80 (ET) 

100 (ET) 


vs. CABR 

- (LTJt) 

20 (LTJt) 

64 (LT.et) 

5 (LTJt) 

33 (LTJt) 

100 (LT.et) 

27 (LTJt) 

41 (LTJt) 

100 (LT.et) 


M20 

25 (LT.et) 

36 (ET) 

-(ET) 

54 (LT.et) 

65 (ET) 

85 (ET) 

60 (LT.et) 

92 (ET) 

100 (ET) 


vs. GINI 

2 (LTJt) 

32 (LTJt) 

100 (LT.et) 

21 (LTJt) 

40 (LTJt) 

100 (LT.et) 

30 (LTJt) 

42 (LTJt) 

100 (LT.et) 




magl.cnt < 21.0 

mag2.cut < 23.0 

mag3.cut < 24.5 


Diagram 

cont. ET 

cont. LT.et 

cont. LT.lt 

cont. ET 

cont. LT.et 

cont. LTJt 

cont. ET 

cont. LT.et 

cont. LTJt 



[%] 

[%] 

[%] 

[%] 

[%] 

[%] 

[%] 

[%] 

[%] 


CCON 

31 (LT.et) 

29 (ET) 

-(ET) 

45 (LT.et) 

41 (ET) 

11 (ET) 

51 (LT.et) 

42 (ET) 

14 (ET) 


vs. CABR 

- (LTJt) 

15 (LTJt) 

81 (LT.et) 

4 (LTJt) 

19 (LTJt) 

86 (LT.et) 

18 (LTJt) 

22 (LTJt) 

98 (LT.et) 


CCON 

29 (LT.et) 

40 (ET) 

2(ET) 

34 (LT.et) 

41 (ET) 

8 (ET) 

48 (LT.et) 

55 (ET) 

28 (ET) 

S 

vs. GINI 

0.5 (LTJt) 

15 (LTJt) 

83 (LT.et) 

2 (LTJt) 

22 (LTJt) 

93 (LT.et) 

6 (LTJt) 

26 (LTJt) 

94 (LT.et) 

X 

CCON 

52 (LT.et) 

69 (ET) 

27 (ET) 

56 (LT.et) 

77 (ET) 

64 (ET) 

71 (LT.et) 

77 (ET) 

73 (ET) 

D 

vs. M20 

4 (LTJt) 

20 (LTJt) 

100 (LT.et) 

11 (LTJt) 

23 (LTJt) 

100 (LT.et) 

13 (LTJt) 

26 (LTJt) 

100 (LT.et) 

S 

GINI 

24 (LT.et) 

36 (ET) 

-(ET) 

30 (LT.et) 

33 (ET) 

6(ET) 

47 (LT.et) 

42 (ET) 

13 (ET) 


vs. CABR 

- (LTJt) 

18 (LTJt) 

88 (LT.et) 

2 (LTJt) 

21 (LTJt) 

94 (LT.et) 

11 (LTJt) 

25 (LTJt) 

100 (LT.et) 


M20 

38 (LT.et) 

46 (ET) 

-(ET) 

48 (LT.et) 

56 (ET) 

34 (ET) 

56 (LT.et) 

49 (ET) 

28 (ET) 


vs. CABR 

- (LTJt) 

18 (LTJt) 

98 (LT.et) 

10 (LTJt) 

25 (LTJt) 

100 (LT.et) 

47 (LTJt) 

24 (LTJt) 

90 (LT.et) 


M20 

35 (LT.et) 

55 (ET) 

10 (ET) 

38 (LT.et) 

63 (ET) 

31 (ET) 

60 (LT.et) 

67 (ET) 

53 (ET) 


vs. GINI 

2 (LTJt) 

27 (LTJt) 

100 (LT.et) 

6 (LTJt) 

24 (LTJt) 

100 (LT.et) 

10 (LTJt) 

30 (LTJt) 

100 (LT.et) 




magl.cnt < 21.0 

mag2.cut < 23.0 

mag3.cut < 24.0 


Diagram 

cont. ET 

cont. LT.et 

cont. LT.lt 

cont. ET 

cont. LT.et 

cont. LTJt 

cont. ET 

cont. LT.et 

cont. LTJt 



[%] 

[%] 

[%] 

[%] 

[%] 

[%] 

[%] 

[%] 

[%] 


CCON 

- (LT.et) 

-(ET) 

-(ET) 

- (LT.et) 

2(ET) 

-(ET) 

6 (LT.et) 

17 (ET) 

-(ET) 


vs. CABR 

- (LTJt) 

7 (LTJt) 

36 (LT.et) 

- (LTJt) 

9 (LTJt) 

31 (LT.et) 

- (LTJt) 

15 (LTJt) 

62 (LT.et) 


CCON 

- (LT.et) 

-(ET) 

-(ET) 

4 (LT.et) 

11 (ET) 

-(ET) 

11 (LT.et) 

30 (ET) 

-(ET) 


vs. GINI 

- (LTJt) 

11 (LTJt) 

48 (LT.et) 

- (LTJt) 

12 (LTJt) 

40 (LT.et) 

- (LTJt) 

19 (LTJt) 

75 (LT.et) 

C 

CCON 

- (LT.et) 

-(ET) 

-(ET) 

7 (LT.et) 

15 (ET) 

-(ET) 

22 (LT.et) 

38 (ET) 

-(ET) 

0 

vs. M20 

- (LTJt) 

3 (LTJt) 

17 (LT.et) 

- (LTJt) 

6 (LTJt) 

20 (LT.et) 

- (LTJt) 

12 (LTJt) 

47 (LT.et) 

s 

GINI 

- (LT.et) 

-(ET) 

0(ET) 

2 (LT.et) 

4(ET) 

-(ET) 

13 (LT.et) 

30 (ET) 

-(ET) 

M 

vs. CABR 

- (LTJt) 

28 (LTJt) 

100 (LT.et) 

- (LTJt) 

35 (LTJt) 

100 (LT.et) 

- (LTJt) 

34 (LTJt) 

100 (LT.et) 

O 

M20 

- (LT.et) 

-(ET) 

-(ET) 

2 (LT.et) 

5 (ET) 

-(ET) 

15 (LT.et) 

32 (ET) 

-(ET) 

s 

vs. CABR 

- (LTJt) 

8 (LTJt) 

30 (LT.et) 

- (LTJt) 

12 (LTJt) 

47 (LT.et) 

- (LTJt) 

12 (LTJt) 

77 (LT.et) 


M20 

5 (LT.et) 

7 (ET) 

-(ET) 

24 (LT.et) 

55 (ET) 

-(ET) 

35 (LT.et) 

75 (ET) 

-(ET) 


vs. GINI 

- (LTJt) 

10.5 (LTJt) 

43 (LT.et) 

- (LTJt) 

17 (LTJt) 

63 (LT.et) 

7 (LTJt) 

26 (LTJt) 

100 (LT.et) 


logCABR 

- (LT.et) 

-(ET) 

-(ET) 

1.3 (LT.et) 

3 (ET) 

-(ET) 

13 (LT.et) 

26 (ET) 

-(ET) 


vs. logASYM 

- (LTJt) 

13 (LTJt) 

45 (LT.et) 

- (LTJt) 

18 (LTJt) 

55 (LT.et) 

- (LTJt) 

24 (LTJt) 

80 (LT.et) 


logCCON 

- (LT.et) 

-(ET) 

-(ET) 

5.5 (LT.et) 

11.3 (ET) 

-(ET) 

11 (LT.et) 

17 (ET) 

-(ET) 


vs. logSMOOTH 

- (LTJt) 

6.3 (LTJt) 

38 (LT.et) 

- (LTJt) 

10 (LTJt) 

34 (LT.et) 

- (LTJt) 

16 (LTJt) 

60 (LT.et) 


logSMOOTH 

56 (LT.et) 

28 (ET) 

64 (ET) 

35 (LT.et) 

39 (ET) 

-(ET) 

42 (LT.et) 

47 (ET) 

46 (ET) 


vs. logASYM 

- (LTJt) 

11 (LTJt) 

- (LT.et) 

- (LTJt) 

18 (LTJt) 

78 (LT.et) 

5 (LTJt) 

22 (LTJt) 

100 (LT.et) 


ASYM 

7 (LT.et) 

6(ET) 

-(ET) 

16 (LT.et) 

31 (ET) 

-(ET) 

28 (LT.et) 

60 (ET) 

40 (ET) 


vs. GINI 

- (LTJt) 

20 (LTJt) 

67 (LT.et) 

- (LTJt) 

25 (LTJt) 

84 (LT.et) 

4 (LTJt) 

30 (LTJt) 

100 (LT.et) 


* Column description: Diagram - commonly used morphological diagrams (for the description of each parameter see Sec. ED; mag.cut - for each survey, 
and in each diagram, the contamination levels are provided at three corresponding magnitude cuts (see Sec. 13.3L cont. ET, cont. LT.et, and cont. LTJt - 
contamination level of ET, LT_et, and LTJt simulated local galaxies, respectively, with the other two morphological types at the corresponding magnitude 
cuts. 
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~ 15% of misclassified sources, and was used in other works, as 
mentioned above. From the established limits we can see again 
how sensitive the diagram is to spatial resolution and data depth, 
changing significantly between the ground- and space-based 
surveys. While in COSMOS the contamination by normal galaxies 
in the region assigned for major mergers, does not exceed the 
values measured by the authors (even at higher redshifts), in the 
analysed ground-based surveys the contamination can be > 70%. 

When applying the M20 - GlNl diagram to classify galaxies 
and select mergers, we highly recommend to use it always with 
additional criteria when dealing with ground-based data, among 
which: a combination with other morphological diagrams, a care¬ 
ful selection of the control sample and evaluation of the diagram 
effectiveness, and/or application of probabilistic approaches. In 
space-based surveys like COSMOS to distinguish between the ET 
and LT sources additional criteria should be taken into account at 
least for galaxies fainter than F814 = 23.0. 


Finally, different combinations of 2 - 3 parameters analysed 
in the previous sections, were used in many studies to distinguish 
between the ET and LT galaxies, and also to study the structure 


nuclear activity, etc (e.g. Pierce et alj 2 

107; Cassata et al. 

2007 

[Bussmann et alj 20091; Povic et al. 2009al 

bl Pentericci et al 

201 c 


IChung et alj|2013l : iBdhm et al.ll2014ll . The combination between 
the concentration parameters (CABR and/or CCO N) and ASYM 


was a ls o use d for selecting ET and L T sources ([ Abraham et al 


1^94, Il996l: I Conselic e et al.l l2006l : iHuertas-Companv et al 
2007 : IPovic et ^ 2009 ah . while in some cases the combina¬ 


tion between the GlNl and CCON w as used to study th e host 
properties of AGN and quasars (e.g. lUrrutia et alJ l2008h . The 
CCON - ASYM - SMOOTH diagram was used i n classifying 
both galaxies and mergers dConselice et alJ l2003l : ICassata et al.l 


th p _^ 

l200^l(^nselice et al.ll200}j) . We found SMOOTH to be especially 
sensitive to both spatial resolution and depth, and the most unstable 
parameter of all analysed. Moreover, together with ASYM, we 
also found it to be more sensitive to noise, as showed in Sec. 13.31 
As already mentioned, we do not recommend the use of these two 
parameters in the ground-based surveys if the noise effects are not 
taken into account carefully. In surveys similar to COSMOS, they 
can be used in combination with other parameters but for sources 
brighter than F814 = 24.0. The combination of only ASYM and 
SMOOTH gives high contamination levels even in the case of 
COSMOS (as showed in Fig.[T3and Table|3- 

We suggest the consultation of Table |2] and Eig. [TT] to 
[T7] when using diagnostic diagrams in morphological studies, 
especially if they are based on 2-3 parameters. We discussed in 
Sec. 14.31 the effectiveness of a particular diagnostic diagram in 
the ground-based surveys and noted that in cases similar to both 
ALHAMBRA and SXDS, special care has to be taken into account 
when classifying galaxies at all magnitude cuts (redshifts). Our 
suggestion when using ground-based data is to avoid the diagnostic 
diagrams based on 2-3 parameters, and to apply instead all 
morphological parameters simultaneously, and to use the statistical 
approaches providing a probability for a galaxy to belong to a 
given morphological type (e.g., approaches similar to that used 
in the galSVM code and tested in Huertas-Company et al. 2008). 
Diagrams that include M20 moment of light showed to be more 
affected by spatial resolution, and less effective in separating 
galaxies. Moreover, the noise should be taken into account if using 
ASYM and SMOOTH parameters. 

In space-based surveys similar to COSMOS, even at the 



Figure 18. Histograms of the di (left panel) and A 2 (right panel) eigen¬ 
vectors resulting from the LDA of the local sample, without including the 
SMOOTH and ASYM parameters. The red cross-hatched, blue hatched, 
and green solid histograms con'espond to the ET, LT_et, and LTJt galaxies, 
respectively. As it can be seen, A[ is the main morphology discriminator, 
whereas A 2 only slightly improves the classification for the LT_et galaxies. 


highest magnitude cuts (F814<24.0), in all diagnostic diagrams 
the contamination levels are significantly lower than at the bright¬ 
est magnitudes in the ground-based surveys. The observational 
bias is insignificant for all parameters at magnitudes brighter than 
F814 = 21.0 (except for the youngest spiral and irregular galaxies), 
as shown in Sec l4.3l however when going to fainter magnitudes the 
contamination increases. Therefore, even when dealing with data 
sets similar to COSMOS, we remind that eventual contamination 
levels are present at fainter magnitudes; the diagnostic diagrams 
that show to be more effective in this case are CCON - CABR, 
CCON-GlNl, CCON-SMOOTH, and/or CABR-ASYM. To 
apply again, when possible, different diagnostics will minimise the 
effect of the observational bias. 


5.4 Results from Linear Discriminant Analysis 

In order to merge all the analysed morphological parameters, and 
to provide an additional and statistical test of our previous results, a 
Linear Discriminant Analysis (LDA; as implemented by Pedregosa 
et al. 2011) was performed. The LDA method allows us to find 
the linear combination of measured parameters that best separates 
two or more classes of objects or events. We tested the combina¬ 
tion of the four morphological parameters (CABR, GINI, CCON, 
and M20) to distinguish between ET, LT_et, and LTJt galaxies in 
the three analysed non-local surveys (ALHAMBRA, SXDS, and 
COSMOS). Knowing that ASYM and SMOOTH are found to be 
affected by noise, we performed the LDA without these two pa¬ 
rameters. However we also performed a LDA including all the six 
parameters for comparison and to check the effect of noise since 
both ASYM and SMOOTH parameters resulted to be the ones more 
significantly affected by noise, mainly in the two ground-based sur¬ 
veys, as we described in Sec. 13.31 

The training phase of the LDA was obtained by fitting the co¬ 
efficients over the measured parameters of the galaxies of the local 
sample. We found that all the possible discrimination is based on 
two eigenvectors (Ji and ^ 2 ) both being linear combinations of the 
four parameters entering in the classification. The first one, di, is 
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Table 3. LDA without SMOOTH and ASYM. 


ALHAMBRA 

Ngal 

score 

SXDS 

Ngal 

score 

20.0 

139470 

0.583 

21.0 

13872 

0.478 

21.5 

139046 

0.510 

23.0 

14579 

0.533 

23.0 

137463 

0.356 

24.5 

14427 

0.525 


COSMOS 

Ngal 

score 

21.0 

49890 

0.549 

23.0 

50150 

0.601 

24.0 

48519 

0.505 


dominated by GINI, CCON, and CABR, whereas the second one, 
A 2 , is based on CABR, GINI, and M20. Almost all the power of 
the discrimination arise from the Ai eigenvector whereas A 2 only 
slightly helps to the classification of LT_et galaxies, as it can be 
seen in Fig. [TS] where the distributions of both di and A 2 for the 
three morphological types are shown. The three concentration in¬ 
dexes show to be more stable, discriminating better between dif¬ 
ferent morphological types in comparison with M20. An attempt 
to improve the LDA was also carried out by doing a second LDA 
iteration including only those galaxies that show the same input 
and predicted output morphology. We did not find any significant 
change, and therefore did not proceed further with this additional 
step. As a result of the LDA of the local sample we found that about 
30% of the objects show a different input morphology than the one 
predicted by the resulting eigenvectors of the LDA. This fact could 
be caused by errors in the measured parameters or even in the orig¬ 
inal classification. To illustrate the results obtained by the LDA we 
present in Figure [TSl the histogram of the predicted morphologies 
after applying the eigenvector discrimination over the local sam¬ 
ple, where for each output morphological type (FT, LT_et, LTJt) 
we show the contribution of galaxies with different input morphol¬ 
ogy. As it clearly can be seen, whereas FT galaxies tend to be well 
classified, LT_FT spread between adjacent input types, and LTJ_T 
are split between both input late-types. 

After obtaining the eigenvectors that better discriminate the 
morphology in the local sample, LDA transformations to the (Ti 
and A 2 ) space were performed for each non local sample (ALHAM¬ 
BRA, SXDX, COSMOS). The results so obtained at the three mag¬ 
nitude cuts and for each non-local sample are presented in Table[3] 
In this table, for each survey, the first column represents the three 
corresponding magnitude cuts, the second column corresponds to 
the total number of galaxies (Ngal) at each magnitude cut and sam¬ 
ple (i.e. galaxies with valid measured parameters after applying the 
filtering criteria explained in Sec. l4.3t while the score column rep¬ 
resents the ratio between the galaxies having the same LDA out¬ 
put morphology as the input one normalized to the total number 
of galaxies Ngal; a score of 1 would indicate a perfect LDA clas¬ 
sification. We also tested the results when including ASYM and 
SMOOTH. In table|4]we present the results obtained from the LDA 
taking into account all morphological parameters. In general, the 
scores obtained after removing these two parameters (see Table 
O are found to be similar to those that would be obtained includ¬ 
ing also ASYM and SMOOTH, however the total number of valid 
galaxies increase significantly and this is more notable in the case 
of the ground base surveys (ALHAMBRA and SXDS). If ASYM 
and SMOOTH are excluded, the total number of galaxies with good 
morphology is found to be larger, thus the final classification get 
worst when including these two parameters as they are indeed af¬ 
fected by noise. Only in surveys similar to COSMOS, the inclusion 
of ASYM and SMOOTH would have a slightly positive net effect. 

Table [ 5 ] corresponds to the obtained scores (ratio of galaxies 
with same output/input morphology over total galaxies) for each 



ET LT_et LTJt 

Predicted Morphology 


Figure 19. Histogram of the predicted morphology obtained from the LDA. 
For each predicted moiphology, the number of galaxies of each different 
input morphological type are displayed. Colours and representations are 
the same as in Fig. d ET galaxies are found to be well classified, however 
LT_ET spread between adjacent input types, and LTXT galaxies ai'e split 
between both late-types. 


Table 4. LDA using all parameters. 


ALHAMBRA 

Ngal 

score 

SXDS 

Ngal 

score 

20.0 

114271 

0.539 

21.0 

10752 

0.491 

21.5 

88186 

0.525 

23.0 

11067 

0.512 

23.0 

66978 

0.479 

24.5 

10077 

0.489 


COSMOS 

Ngal 

score 

21.0 

49292 

0.653 

23.0 

47988 

0.655 

24.0 

43061 

0.600 


morphological type and magnitude cut in each survey, as resulting 
from the LDA with neither SMOOTH nor ASYM included. As can 
be seen, in all surveys, FT galaxies are better recognised at the 
brighter magnitude cuts, whereas LTs tend to be better recognised 
at fainter magnitude cuts. 


6 SUMMARY AND CONCLUSIONS 

In this paper we present for the first time a systematic study of 
the impact from spatial resolution and depth on the six morpho¬ 
logical parameters commonly used for galaxy morphological 
classification: Abraham concentration index, Gini coefficient, 
Conselice-Bershady concentration index, M20 moment of light, 
asymmetry, and smoothness. We studied how strong the impact 
is for three data sets that have different spatial resolutions and 
depth: ALHAMBRA and SXDS, as an example of ground-based 
data, and COSMOS, as an example of deep space-based data. Our 
results correspond to maximum analysed redshift values of: 1.0, 
2.2, and 1.6 in the three surveys, respectively, where the provided 
values describe the distribution of 95% of each sample. We used 


Table 5. LDA results by morphological type. 


ALHAMBRA 

ET 

LT.ET 

LT.LT 

SXDS 

ET 

LT.ET 

LT.LT 

20.0 

0.86 

0.51 

0.02 

21.0 

0.73 

0.33 

0.15 

21.5 

0.47 

0.74 

O.IO 

23.0 

0.72 

0.51 

0.10 

23.0 

0.17 

0.55 

0.41 

24.5 

0.55 

0.67 

0.13 




COSMOS 

ET LT.ET 

LT.LT 






21.0 

0.99 

0.29 

0.00 






23.0 

0.83 

0.54 

0.15 






24.0 

0.55 

0.55 

0.15 
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a sample of 3000 early- and late-type local galaxies, with the 
available visual morphological classification, and we measured 
their parameters in 2 cases: 

- first, the reference values, that correspond to the real galaxy 
redshifts (magnitudes), and 

- second, the simulated values, obtained after moving the local 
galaxies to the observational conditions and magnitude/redshift 
distribution of three selected surveys. 

Comparing the results obtained in these two cases we showed 
how each parameter changes in each survey, at three particular 
magnitude cuts. Finally, we analysed and quantified how the 
impact from spatial resolution and depth alfects some of the most 
used morphological diagrams, gave some suggestions for the 
galaxy/merger classification studies, and used the LDA statistical 
approach to analyse the most effective combination of parameters 
to distinguish between different types. 

In the following, we summarise some of our main findings: 

- All six analysed morphological parameters suffer from sig¬ 
nificant biases related to the spatial resolution and data depth. 

- The impact of the spatial resolution on the morphology is 
much stronger in comparison with the data depth, being therefore 
the most responsible for changing the parameters in the ground- 
based surveys, making in general the galaxies to appear less con¬ 
centrated and more symmetric. 

- We stress that ASYM and SMOOTH are more sensitive to 
noise effects. Survey noise should be taken into account carefully 
when using these two parameters in morphological classification of 
galaxies, especially when dealing with ground-based data. 

- M20 results to be also significantly affected in all surveys, 
changing both the shape and range of its distribution at all bright¬ 
ness levels. However, M20 - GINI is the most used diagram for se¬ 
lecting interacting (merging) systems. We highly recommend the 
use of this diagram simultaneously with other morphological pa¬ 
rameters, when dealing with ground-based data sets, and in space- 
based surveys like COSMOS at least when studying sources fainter 
than F814 = 23.0. 

- CCON shows to be more sensitive to the spatial resolution 
in comparison to GINI and CABR. It works therefore better for 
space-based data, whereas the other two concentration indexes be¬ 
have better for ground-based surveys. 

- In surveys similar to ALHAMBRA, when analysing the 
highest density regions of early- and late-type galaxies in the main 
morphological diagnostic diagrams, the impact from spatial resolu¬ 
tion and data depth introduces contamination levels of 5 - 25% for 
ET and 3 - 35% for LT galaxies (depending on the diagram), for 
the brightest galaxies with F613 ^ 20.0. The diagrams that seem to 
work the best in this case are those that do not include M20 moment 
of light. At the faintest analysed magnitudes (F613 ^ 23.0) the con¬ 
tamination levels increase significantly, being as high as 60 -100%, 
making each of the diagnostics useless if used separately. Similar 
results are obtained in the case of SXDS, but at higher magnitude 
cuts (of order 1 - 2). Moreover, in both surveys, at all magnitude 
cuts, classification ot LTJt galaxies suffers very high levels of con¬ 
tamination. Taking all this into account, when dealing with ground- 
based data sets, we suggest to avoid the use of 2 - 3 parameter di¬ 
agnostic diagrams in morphological classification, and to apply in¬ 
stead the use of all morphological parameters simultaneously, and 
to use statistical approaches based on probability distributions that 
the galaxy is ET or LT. 

- In space-based surveys similar to COSMOS, even at the 
highest magnitude cuts (F814 < 24.0), the contamination levels are 


significantly lower than those at the brightest magnitudes in the 
ground-based surveys. The observational bias is insignificant for 
all parameters at magnitudes brighter than E814 = 21.0 (except for 
the LTJt galaxies). However, when going to fainter magnitudes the 
contamination increases, and depending on the diagram it changes 
from 2-25% at F814 < 23.0, to 6 - 35% for ET, and 20 - 70% for LT 
galaxies at F814 < 24.0. In surveys similar to COSMOS, we again 
suggest to use several diagnostics to classify galaxies when going 
to magnitudes fainter than F814 = 23.0. 

- Through LDA analysis we obtained that the combination of 
CABR, GINI, and CCON parameters, is the most effective to dis¬ 
tinguish between different morphological types. 

The results presented in this paper can be directly applied to 
any survey similar to ALHAMBRA, SXDS and COSMOS, and 
also can serve as an upper/lower limit to take into account when 
classifying galaxies using shallower/deeper data sets. 
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